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ABSTRACT 

A number of ultra-cool dwarfs emit circularly polarised radio waves generated by the electron 
cyclotron maser instability. In the solar system such radio is emitted from regions of strong auroral 
magnetic field-aligned currents. We thus apply ideas developed for Jupiter's magnetosphere, 
being a well-studied rotationally-dominated analogue in our solar system, to the case of fast- 
rotating UCDs. We explain the properties of the radio emission from UCDs by showing that it 
would arise from the electric currents resulting from an angular velocity shear in the fast-rotating 
magnetic field and plasma, i.e. by an extremely powerful analogue of the process which causes 
Jupiter's auroras. Such a velocity gradient indicates that these bodies interact significantly with 
their space environment, resulting in intense auroral emissions. These results strongly suggest 
that auroras occur on bodies outside our solar system. 

Stars: magnetic field, late-type, low-mass - Planets and satellites: 



polarised nature of the emission implicates the 
electron cyclotron m aser instab i lity (CMI) as the 



Subject headings: brown dwarfs 
magnetic fields, aurorae 

Introduction 



Ultra-cool dwarfs (UCDs) are objects with 
spectral type later than M7, comprising the lowest 
mass stars and brown dwarfs. Twelve of these (out 
of ^200 observed) have been found to be intense 
sources of radio emissions with spect ral luminosi- 

(iBerger et al 



ties t y picall y of order a MW Hz 
2001 1 l2009h iBergerl [200l iHallinan et all 12 006. 



2007 ,1 20081; I Antonova et~alll2007il2008HPhan-Bao et al 



2007:lMcLean et al.ll201lLl2012tlRoute fc Wolszczan 



2012). The unpolarised component probably in- 
cludes synchrotron emission, but the radio emis- 
sion from a number of these, which have fast 
(^2 h) rotation periods and strong (^0.1 T) mag- 
netic fields, has been shown to be highly circu- 
larly polarised and modulated at the bod ies' ro- 
tation periods (e.g. Hallinan et al. 20061 ). The 



sourc e mechanism (IWu fc Led Il979t Treumann 
20061) . lYu et all (|201lJ) modelled the radio emis- 
sion, assuming it is generated by unaccelerated 
precipitating hot plasma trapped on coronal loops. 
However, this process would be likely limited to 
short-lived flare-type events, does not directly re- 
sult in increas ed radio power with faster rotation 
as is observed ( McLean et al. 20121 ). and the CMI 
is known to be generated at much greater efficiency 
when electrons are accelerated dow n magnetic 
field lines by field-aligned voltages ([Treumann 
20061 ). The presence of sustained CMI-generated 



radio emissions is thus instead strongly suggestive 
of the existence of quasi-steady auroral magnetic 
field-aligned currents. Indeed, all CMI source 
regions observed in situ at planets in the solar 
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system exhibit accelerated electron populations, 
and radio emissions generated by the CMI have 
been extensively shown to be closely associated 
with auroral emission, caused when downward- 
precipitating ele c trons impact the atmosphere 
(iTreumanij l2006t IClarke et all [20091: lLamv et all 
20091: lLamv et al.1 1201ft iNichols et al.1 l2010aM . 
Yu et all (|2012h have recently developed their 
analysis to consider the results of an accelerated 
population of electrons, and their simulations in- 
dicate that such populations could produce the 
observed UCD emissions. 

Field-aligned currents arise from a divergence 
in field-perpendicular currents, which are, through 
E = —v x B, driven by plasma velocities relative 
to the neutrals in the conducting outer layer of 
the atmosphere. Thus, strong field-aligned cur- 
rents flow when there is a sharp gradient in this 
departure from corotation, giving rise to a strong 
divergence in the field-perpendicular currents. We 
may therefore infer from the observation of CMI- 
generated radio emissions from UCDs that such 
angular velocity gradients with auroral currents 
exists in the magnetospheres of these objects. At 
planets in the solar system angular velocity gradi- 
ents occur near the boundary between open and 
closed field lines, as at Earth and Saturn, or are 
due to centrifugally-driven outfl ow of internally - 
generated plasma, as at Jupiter. Schriiverl ([2009) 
presented the hypothesis that the radio emission 
from UCDs could result from a Jupiter-like cur- 
rent system on the basis of a simple dimensional 
scaling relation, but did not consider quantita- 
tively the magnetosphere-ionosphere (M-I) cou- 
pling current system. Here we have developed a 
simple axisymmetric model of the currents aris- 
ing from departure from rigid corotation in the 
magnetospheres of UCDs, based on tha t used pre- 
vious l y to study Jupiter's au r oral oval ( HilJ 197S 
20011 ICowlev fc Buncd boOll ICowlev et all 12002 



Nichols fc Cowlevll2003l 12004 120051; ICowlev et all 
20051 : iNicholsl l2011bh and to estimate the radio 
luminosity of fast-rot ating Jupiter-like exoplanets 
([Nichols! 1201 la , 2012 ). We explain the properties 
of the radio emission from UCDs by showing that 
it would arise from the electric currents resulting 
from an angular velocity shear in the fast-rotating 
magnetic field and plasma, i.e. by an extremely 
powerful analogue of the process which causes 



Jupiter's or Saturn's auroras. Such a velocity 
gradient indicates that these bodies interact sig- 
nificantly with their space environment, resulting 
in intense auroral emissions. 



2. Theoretical background 
2.1. Basic assumptions 

For the purposes of development of the model, 
we assume that a jovian picture of an ultra-cool 
dwarf (UCD) is appropriate, i.e. that the UCD's 
magnetic field is generated deep within its interior, 
such that the external field generated by the body 
is mainly dipolar, and that a non-conducting at- 
mosphere surrounds the body, the upper portion 
of which is maintained to some degree conductive 
by photon and/or corpuscular impact from the 
outside. The atmospheres of cool UCDs are not 
expected to be strongly ionised, although it has 
been suggested that dust clouds may indu ce some 
degree of ionisation ( Helling et al. 2011al fb). At 
solar system planets a number of sources of con- 
ductivity are known besides solar X-ray and EUV 
radiation, such as galactic cosmic rays and au - 



roral electron precipitation ( Rvcroft et al. 2008) 



The latter is particularly important, since it re- 
sults in a positive feedback which greatly increases 
the intensity of the auroral currents. At Jupiter, 
for example, modeling work indicates that the 
feedback from auroral precipitation increases the 
conductance from a residual value of ^0.01 mho 
(where 1 mho = 1 Siemens) generated e.g. by so- 
lar X-ray and EUV irradiation to values of up to 
10 m ho ([Strobel fc Atreval Il983t iMillward et al 



2002), such that the feedback plays a dominant 
role in influencing the lo cation and intensity of th e 
M-l coupling currents (|Nichols fc Cowlevl 120041) . 
Thus we suggest that given some small initial con- 
ductivity generated by a source mentioned above, 
the auroral current system would amplify itself via 
the feedback from the auroral precipitation. The 
exact details of such feedback requires ionospheric 
modeling that is beyond the scope of the present 
paper, but suffice it to say that due to this effect 
the low temperatures of UCDs do not necessarily 
prohibit the flow of strong M-I coupling currents. 
In the absence of such an ionospheric model for 
UCDs, in the present work we simply take the 
Pedersen conductance £p to be a constant, dis- 
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cussed further below. 

A principal assumption is that the spin and 
magnetic axes are roughly co- aligned, since this 
is the configuration that has been most studied 
for bodies in our solar system, although we note 
that it is probable that only ^50% of UCDs ex- 
hibit this polarity. With reversed field polarity, 
the associated M-I coupling current system de- 
scribed below would also be reversed, such that 
the upward current, which in Fig. 1 is strongly 
peaked at 15°, would instead be widely distributed 
at high latitudes, and the strong downward cur- 
rents would be easily carried by upward-going 
ionospheric electrons, leading to little auroral and 
radio emissions. We also assume approximate ax- 
isymmetry about the magnetic axis, such that 
the flows in the ionosphere are considered to be 
wholly azimuthal about that axis, and can thus be 
described in terms of departure from corotation 
with the neutral atmosphere as a function of colat- 
itude 0i from the magnetic axis. Such a simplifica- 
tion seems appropriate for representing to lowest 
order the dynamics of a fast rotating magneto- 
sphere, although we note that some processes re- 
lated to an open magnetosphere will consequently 
be excluded from the model, such as 'cusp-like' or 
'substorm-like' processes which would be ordered 
in azimuth with respect to the object's velocity 
through the surrounding medium. We further as- 
sume for simplicity that departures of the upper 
neutral atmosphere from rigid corotation due to 
ion-neutral drag also take the form of winds that 
are approximately axisymmetric about the mag- 
netic axis, thus neglecting the effect of the Coriolis 
force on such motions in cases where the spin and 
magnetic axes are significantly inclined. 



It has been suggested by iKuznetsov et al 



(|2012n that the short duration of the spikes in 
the CMI-induced radio emissions is indicative of 
an active sector in the magnetosphere of a UCD. 
We suggest, however, that the observed short 
duty cycles of < 0.1 5 for UCD radio emissions 
(Halli nan et al.ll2006l ) are not inconsistent with an 
essentially axisymmetric field-aligned current dis- 
tribution. Jupiter's main auroral oval is to first 
order described by an axisymmetric annulus of 
isotropic emission ~1° wide located ~15° from the 
magnetic pole, which is tilted by ~10° from the 



spin axis (iGrodent et al]|2003l: Iciarke et al.ll2004 ; 
Nichols et al.l 120091 ) . Therefore, as the dipole axis 
cones around the spin axis, from a viewpoint near 
the planet's spin equatorial plane over half of the 
auroral oval is visible o ver roughly a third of the 
planet's rotation period (|Nichols et al. I l2007l) . The 
associated radio emission, however, is strongly 
beamed, into ~ 1 .6 sr in the case of the HOM and 
DAM emissions ( Zar ka et al . 2004). Hence, the 
radio emission peaks strongly as the beam rotates 
into view, and the FWHM of the occurrence prob- 
ability distribution of the 'Source A' DAM emis- 
sion (i.e. that used to determine Jupiter's System- 
Ill rotation period of ~9 h 55 min ) is ~50°, i.e. 
an a pparent duty cycle of ~0.14 (Hig gins et al.l 
1996(). This is consistent with the duty cycles of 
< 0.15 for UCDs. Having said this, it is important 
to recognise that axisymmctry is simply a lowest 
order approximation. Saturn's radio emissions, 
for example, are observed from all loc al times but 



are b rightest in the morning sector (jLamv et al 



2009), despite the internal planetary magnetic 
field being axisymm etric to within detection limits 
(Bur ton et al. Il2010h. F i nally, as has also been sug- 
gested by Berger et al. ( 20091) . it is probable that 
as is observed at Jupiter the magnetic poles may 
not lie on respective antipodes, due to the non- 
dipolar terms in the internal field, such that the 
radio emission from the two hemispheres would 
not necessarily be expected to be observed 180° 
apart in phase. 



2.2. Current system equations 

We consider three angular velocities in the 
computation of the M-I coupling currents. The 
first is the angular velocity of the deep inte- 
rior of the UCD ilucD, taken initially to be 
8.7 x 10 -4 rad s _1 in conformity with the observed 
periods of ^2 h for radio-bright UCDs. The sec- 
ond is the angular velocity of the plasma above the 
ionosphere uj, taken to be constant on each flux 
tube, and the third is the angular velocity of the 
atmospheric neutrals in the Pedersen layer of the 
ionosphere i7^ C£) , which is expected to lie some- 
where between ftjjCD and ui due to the frictional 
drag of ion-neutral c ollisions ( Huang fc Hill 1989; 
Millward et~al~ll2005h . In such a case we have 
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{Qucd — ^ucd) = k(flucD — w) 



(1) 



where < k < 1. The actual value of k is 
somewhat uncertain for Jupiter, with modelling 
work indicating a value of ~0.5 ([Millward et al. 
2005), which is employed in the previous jovian 



modelling work discussed and also now adopted 
here, although the results are not expected to 
be sensitively dependent on reasonable choices. 
We should note that in the following work we 
consider only sub-rotation of the plasma with re- 
spect to the UCD. For Jupiter, super-rotation 
of the magnetospheric plasma may occur follow- 
ing strong solar wind-induced compressions of the 
magnetosphere, with a correspondi ng reversal of 



the current system descr i bed below (Hanlon et al.. 
2004 ICowlev et alj|2007t iNichols et al.ll2007t) . but 



in the steady state the plasma sub-rotates, associ- 
ated with the transfer of angular momentum from 
the atmosphere to the external medium. 

The equatorward-directed height integrated 
ionospheric Pedersen current is then given by 



ip = T,pEi = T,ppi(fl 



UCD 



w)Bi 



(2) 



where Ei is the electric field in the rest frame of the 
neutral atmosphere, pi — Rjjcd sin#i is the dis- 
tance from the spin axis (where we assume initially 
the UCD has a radius Rucd equal to Jupiter's 
equatorial 1-bar level radius Rj = 71, 373 km since 
all bodies of roughly solar composition from giant 
planets through to the very lowest mass stars have 
radii roughly similar to Jupiter). Employing equa- 
tionQ]in equation [5] and introducing the 'effective' 
Pedersen conductance £p, reduced from the true 
value by a factor of (1 — k) owing to the slippage 
of the neutral atmosphere from rigid corotation, 
yields 



ip = T,* P pi(QucD - oS)Bi 



(3) 



As mentioned above, we simply take a constant 
value for £p, initially 0.5 mho, which is similar to 
the values employed in previous Jupiter modelling 
works, and, as shown below, yields radio power 
values in agreement with observations of UCDs. 
However, in the results below we also examine the 



effect of different values of this parameter. The to- 
tal Pedersen current at a given co-latitude is then 
found by integrating in azimuth, such that 



Ip = 2irpiip = 2irY,ppf({l 



UCD 



(4) 



the divergence of which gives the field- aligned cur- 
rent density just above the ionosphere, i.e. for the 
northern hemisphere 



J\\i 



1 



2npf svaOi d9i 



(5) 



2.3. 



Field-aligned acceleration and energy 
flux 



The upward-directed field-aligned current com- 
puted in the previous section can not, in general, 
be carried simply by precipitating magnetospheric 
electrons alone, and must be driven by down- 
ward acceleration of those electrons by a field- 
aligned voltage. The maximum field-aligned cur- 
rent density that can be carried by an unaccclcr- 
ated isotropic Maxwellian population is 



( w th 



1/2 



(6) 



where e, m e , n and Wth are the charge, mass, 
number density and thermal energy of the electron 
source population, respectively, the latter being 
equal to equal to kpT, where kp is Boltzmann's 
constant and T is the temperature. In common 
with previous works we use the values for the high 
latitude hot plasma at Jupiter, i.e. n = 0. 01 cm~ 3 
and Wth = 2-5 keV (|Scudder et al.lll98ll ). which, 
as shown below yield results consistent with ob- 
servations, but we also examine the effects of wide 
ranges of values for these parameters. The corre- 
sponding unaccelerated kinetic energy flux is 



E f o = lenWth 



W th 
2irm r - 



1/2 



(7) 



For field-aligned current densities larger than j\\io, 
a field-aligned voltage is required. For Earth, the 
current- voltage relation derived using the kinetic 
theory of lKnight (Il973h is applicable, but for more 
powerful systems in which either the source elec- 
tron population is very energetic initially, or be- 
comes so following acceleration by voltages com- 
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parable to or exceeding the electron rest mass en- 
ergy (^511 keV), Cowley's (2006) relativistic ex- 
tension to Knight's theory becomes appropriate. 
In the results shown in this paper it is clear that 
the relativistic theory is required in the case of 
M-I coupling current systems at UCDs, such that 
we employ the relativistic current-voltage relation 
given by 



j\\i< 



= 1 



( m e c 2 \ 
Wth ) 



l 



(8) 



where c is the speed of light and $ is the mini- 
mum voltage required to drive the current jiu at 
ionospheric altitude, applicable in the limit that 
the accelerator is located at high altitude along 
the field line (where B <C Bi). The corresponding 
precipitating electron energy flux is given by 



E f 



E 



fo 



e<E> 



(w th J 


6 


2 


O / m e c 2 


}+3 



(9) 

Equations [6][9j along with values for the electron 
source population number density and tempera- 
ture as discussed further below, are thus used to 
determine the precipitating energy flux resulting 
from the current system discussed above. The to- 
tal precipitating power for each hemisphere P e is 
then obtained by integration of Ef over the hemi- 
sphere, i.e. 



90 

/', = / 2irRl CD sm6 l E f d9 % , (10) 



and we finally obtain the spectral luminosity using 



L f — 



Pe 



100 Av 



(11) 



where we assume that the beam from only one 
hemisphere is observable at any one time and we 
take the electron cyclotron maser instability to 
have a generation efficiency of ~1%, a value ob- 
served at both Jupiter and Saturn. Specifically, 
at Saturn this value was directly measured dur- 
ing the traversal of the Cassini spacecraft through 



the source reg ion of the Saturn Kilometric Radi- 
ation (SKR) (jLamv et al.l 1201 If ) . and at Jupiter 
the radio power output is ~1% of the total pre- 
cipitating electron power, as measured from ob- 
servations of the UV aurora and computed the- 



oretically (iGustin et al.l 12004 IClarke et al.l l2009t 
Cowlev et al.l 2005 ). The bandwidth Av is as- 
sumed to be equal to the electron cyclotron fre- 
quency in the polar ionosphere, i.e. 



eB t 
2irm e 



(12) 



an approximation validated by obse rvation s of so - 
lar system planets. For example, IZarka (Il998h 
argue that the jovian b-KOM, HOM and DAM 
emissions are a single radio component, thus with 
a frequency range of ~10 kHz to ^40 MHz, and 
at Saturn the SKR extends from a few kHz to 
^1200 kHz. This arises because the CMI occurs 
at all points below the field-aligned voltage drop, 
situated at least a few planetary radii up the field 
line, and the top of the ionosphere at 1 plane- 
tary radius. The emission is generated at the local 
electron cyclotron frequency, such that the band- 
width is then determined by the difference between 
the field strengths at these two locations. The 
strength of a dipole field drops off approximately 
as cube of the distance along the field line, such 
that at altitudes of only ~1 planetary radius up 
the field line the frequency drops by an order of 
magnitude. The bandwidth is thus reasonably well 
represented by the high frequency cutoff, which is 
also the frequency used to measure the magnetic 
field strength in the polar region. 

2.4. Angular velocity profile 

We now discuss the model for the plasma an- 
gular velo city used in this pap er. Following previ- 



apc 

ous work (|Cowlev et al.ll20051) , the gradient in the 



angular velocity is conveniently represented by a 
tanh function with respect to co-latitude, given by 



UCD 



UCD J L 

-) ) (l + tanh 



where (uj /Qjjcd)l is the low value to which the 
angular velocity transitions near the magnetic pole 
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Fig. 1. — Figure showing profiles of model M- 
I coupling current system parameters. Specif- 
ically we show (a) the plasma angular velocity 
normalised to that of the UCD (uj/CIucd): (b) 
the azimuth-integrated ionospheric Pedersen cur- 
rent Ip in GA, (c) the field-aligned current den- 
sity at the top of the ionosphere j\u in mA m~ 2 , 

(d) the minimum field-aligned voltage required to 
drive the current plotted in panel (c) in MV, and 

(e) the resulting precipitating electron energy flux 
Ef in kW m -2 , all plotted versus ionospheric co- 
latitude 6i. The Pedersen conductance employed 
is Ep = 0.5 mho and the source electron num- 
ber density and thermal energy are jovian values 
of 0.01 cm -3 and 2.5 keV. The horizontal dot- 
ted lines in panels (a) and (c) indicate values of 
(us / Clue d) — 1 and j\u = mA m -2 , respectively. 



from rigid corotation at lower latitudes, and 9i C 
and A9i are the latitudinal centre and half-width 
of the transition region, respectively. We initially 
employ values for these parameters which give 
a Jupiter-like angular velocity profile, which on 
that planet generates the field-aligned currents 
which cause the main auroral oval and the HOM, 
b-KOM and non-Io-DAM radio emissions, i.e. 
(uj/CI p ) l = 0.25, 6ic = 15° and A0* = 0.5°. The 
actual form of the angular velocity gradient de- 
pends on the physical processes causing the shear, 
as is examined further in the Discussion, and in 
the results below we also examine the effects on 
the radio power of taking different values for these 
parameters. 



3. Results 

3.1. Results with representative values of 
system parameters 

In this section we present results of the model 
using appropriate spot values for the various 
model parameters as discussed above, before go- 
ing on in the next section to examine ranges in 
the parameter space. Figure [T] shows represen- 
tative results, in which we employ a uniform 
ionospheric field strength of 0.3 T and rotation 
period of 2 h, in conformity with typical values 
obtained from the bandwid th and modulation o f 
the UCD radio emissions ( Hallinan et al.l [2008). 
First, Figure [T^, shows the angular velocity pro- 
file employed, given by equation 1131 Moving from 
large to small co-latitudes (or equivalently on field 
lines which thread the equatorial plane at increas- 
ing radial distances), the plasma angular veloc- 
ity initially near-rigidly corotates, and transitions 
over ~1° at 15° co-latitude to a quarter of rigid 
corotation at smaller co-latitudes. The resulting 
equatorward azimuth-integrated ionospheric field- 
perpendicular (Pedersen) current given by equa- 
tion [5] is shown in Figure [TJd, increasing from small 
values near the pole to a value of ^58 GA, before 
falling rapidly due to the plasma angular veloc- 
ity gradient shown in Figure UK- The resulting 
field-aligned current given by equation [5] is shown 
in Figure [TJ;, which is negative, i.e. downward, 
near the pole and switches to a significant peak of 
positive, i.e. upward, field-aligned current values, 
reaching amplitude ^0.7 mA m -2 , centred on 15°. 
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Fig. 2. — Plots showing the spectral luminosity 
computed in the model for different pairs of values 
of the source electron population number density 
n and thermal energy Wth- In panels (a-d) we 
use Pedersen conductances of Y,* p = 0.005, 0.05, 
0.5, and 5 mho, respectively. The crosses indicate 
the jovian values employed in Section [2.3[ and the 
black contour indicates where L r = 1 MW Hz -1 , 
comparable to observed values. 



This is the region of auroral field-aligned current 
that sustains the unstable electron distributions 
responsible for the CMI and thus the radio emis- 
sions. The required field- aligned voltage given by 
equation [8] is then shown in Figure QJI, peaking 
at —11 MV, and the corresponding precipitating 
electron energy flux given by equation [9] is shown 
in Figure QJ, peaking at —9 kW m~ 2 . Integrating 
this energy flux and converting to spectral lumi- 
nosity as discussed above yields a radio power 
output of — 1.1 MW Hz -1 , i.e. consistent with the 
typical observed radio luminosities of UCDs. We 
also note that the bandwidth of the radio emis- 
sion, equal to the electron cyclotron frequency 
of —8.4 GHz, is perforce similar to that observed, 
since the ionospheric field used is determined from 
the measured radio frequencies. 



3.2. Parameter space investigation 

In the above section we have used reasonable 
values for a number of system parameters based 
on observations at Jupiter in order to demon- 
strate the model and indicate how an angular 



velocity gradient leads to radio luminosity that, 
for the case of fast-rotating UCDs, is comparable 
to that observed. However, the only parameters 
appropriate to this model that are known from 
observations, besides the radio power, are the po- 
lar magnetic field strength and rotation rate of 
such objects. The other parameters, such as the 
form of the angular velocity profile, the density 
and temperature of the source electron population 
and the ionospheric conductance are not known. 
Here we therefore examine ranges of these param- 
eters to determine what values are consistent with 
the observed radio emission. We first examine 
the effects of the conductance and source electron 
population parameters, given the angular veloc- 
ity profile given by equation [T3l and we go on to 
investigate the effect of changing the form of the 
angular velocity profile. 

In each panel of Figure [5] we show the radio 
spectral luminosities obtained from the model us- 
ing different pairs of values for the source electron 
thermal energy Wth and number density n. The 
values adopted above obtained from Voyager ob- 
servations of the plasma outside the current sheet 
at Jupiter, i.e. 2.5 keV and —0.01 cm~ 3 , respec- 
tively, are shown by the crosses in each panel. 
The low density arises since the cold plasma is 
centrifugally confined to the equatorial plane, and 
only the hot, rarefied plasma population has a cen- 
trifugal confinement scale hei ght large enough to 
populate the entire flux tube (jHill fe Michellll976t 
Caudal 1986). At high latitudes, the low beta 
plasma, combined with the field-aligned voltage 
and loss of particles to the ionosphere, is then 
favourable for generation of the CMI. We note 
that previous estimates of the electron number 
density and temperature on the flux tubes pro- 
ducing the CMI at UCDs have been obtained, i.e. 
1.25-5 xlO 5 c m" 3 and 1 - 5 xlO 7 K (equivalent to 
-0.86 keV) bv lYu et all (|201ll ). but these assume 
that the CMI is produced by an unaccelerated 
loss-cone population. We suggest, however, that 
at fast rotating UCDs, at which the magneto- 
spheric plasma is likely to be significantly more 
centrifugally confined than at Jupiter, the high 
latitude plasma parameters may possibly be closer 
to the jovian case. In light of the uncertainty sur- 
rounding the source plasma population, we show 
the results of the model using the angular velocity 
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profile given by equation [T3l for wide ranges of 
these parameters, i.e. 5 orders of magnitude ei- 
ther side of the jovian values. In panels (a-d) we 
employ £p = 0.005, 0.05, 0.5, and 5 mho, respec- 
tively, reflecting the uncertainty in the Pcdcrsen 
conductance. The contour labelled "6" indicates 
where the radio luminosity is equal to 1 MW Hz -1 , 
representative of the observed luminosities. 

It is apparent that over the majority of the pa- 
rameter space plotted, the luminosity increases 
with decreasing n and increasing Wth, with a 
change of slope as the relativistic terms dominate 
above Wth — 511 keV. Where the slope reverses in 
the top right corner of each panel the field-aligned 
current can be carried by an unaccelerated pop- 
ulation. For simplicity, in this region we have 
assumed that in the auroral region the loss-cone 
is filled by strong pitch angle diffusion, and that 
the required net current is then created by a cor- 
responding upward flux of electrons. We have also 
taken the same CMI generation efficiency of 1% 
as for the accelerated population, despite the less 
favourable loss-cone distribution in this case. Both 
these assumptions will likely lead to overestimates 
of the radio flux from this region. We also note 
that the previous estimates of the sourc e popula- 
tion parameter values discussed above ( Yu et all 
20111 ) are off the top of each panel, thus clearly 
in the unaccelerated regime. In any case, it is 
apparent in the accelerated region that a wide 
range of parameter values are able to produce 
the required radio luminosity, in general requir- 
ing hotter, less dense plasma for lower Pedcrsen 
conductance. The plot also confirms the results 
shown above, i.e. that typically jovian values for 
the source plasma population parameters and the 
Pedersen conductance lead to observed UCD ra- 
dio luminosities. 

We now consider the effect of the form of the 
angular velocity profile, and show in Figure [3] the 
radio luminosities obtained when the parameters 
in equation [13] are varied individually. Specifi- 
cally, from top to bottom the parameters varied 
are (a) the low value to which the angular veloc- 
ity transitions from rigid corotation (b) 
the latitudinal half-width of the transition region 
A0i, and (c) the centre of the transition region 9i C . 
As a guide, the horizontal dotted lines highlight 




Fig. 3. — Plots showing the spectral luminos- 
ity computed in the model for different values of 
the angular velocity profile parameters. These 
are (a) the low value to which the angular ve- 
locity transitions from rigid corotation (lo/Q. p )l, 
(b) the latitudinal half- width of the transition re- 
gion A0i, and (c) the centre of the transition re- 
gion 8i C . The horizontal dotted lines are located at 
L r = 1 MW Hz -1 . As previously, source electron 
population parameters are fixed at jovian values. 
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Fig. 4. — Figure showing modelled (pluses) 
and observed (crosses) radio luminosities of the 
three confirmed UCD sources of CMI-generated 
radio emissions versus each UCD's angular veloc- 
ity. Model parameters are given in Table [T] 

1 MW Hz -1 , which delimits the parameter ranges 
reasonably consistent with observations. Thus, 
an angular velocity profile which transitions to 
< 70% of rigid corotation across a region of half- 
width < 2.5° at a co-latitude > 9° will produce 
radio powers consistent with observations. These 
are not particularly tight constraints, indicating 
that the model does not require fine tuning to 
produce the required emissions. Future modelling 
of the radio beams, along the lines of previous 
work that has concentrated on localised active 



regions (jKuznetsov et al.ll2012l) . from angular ve- 
locity transition regions of different geometries 
would be useful to constrain the angular velocity 
profiles further. 



3.3. Comparison with observations 

We now compare the model spectral luminosi- 
ties with the confirmed sou rces of CMI-induced ra - 
dio emissions discussed by lHallinan et al. I fl2008h . 
i.e. TVLM 513-46546, LSR J1835+3259, and 
2MASS J00361617+1821104 (hereafter TVLM 
513, LSR J1835 and 2M J0036), which have 
relevant properties listed in Table [1] We esti- 
mate the spectral luminosities of the CMI-induced 
radio emission fro m the peak amplitudes an d 
widths plotted by lHallinan et~all (|2007l [2008). 



As discussed above, CMI-induced radio emission 
is strongly beamed, into 1.6 sr for the case of 
Jupiter's HOM and DAM emissions, resulting in 
a duty cycle of ~0.14 for these emissions. We esti- 
mate the beam width a for these UCDs by simply 
scaling to the jovian emissions via 



Ar/r 
0.14 



1.6 



(14) 



where the pulse duty cycle Ar/r is given in Ta- 
ble [1] This, however, does assume that the angular 
velocity profile, and thus that of the field-aligned 
currents, are the same for each dwarf, which may 
not be the case. The spectral luminosity then fol- 
lows from the peak flux density F via 



L r = Fad 2 



(15) 



where d is the distance to the object from the 
Earth. This is shown by the crosses in Figure |U 
In comparison, the modelled spectral luminosity, 
shown by the pluses, is computed using the radii, 
rotation periods, polar ionospheric field strengths 
and Pedersen conductances as listed in Table [T] 
The Pedersen conductances employed for TVLM 
513 and 2M J0036 are 0.5 mho, as is also used 
in Figure [TJ but it was found that 0.3 mho yields 
better agreement for LSR J1835. It is also worth 
noting that employing 0.8 mho for all three dwarfs 
yields spectral luminosities consistent to within 
~ 7% of the quiescent valu es (computed assumin g 
isotropic emission) given bv lHallinan et al.l (]2008f h 



We also note the potential CMI-induced emis- 
sions fro m the L dwarf binar y 2MASSW J0746425+ 
200032 (jBerger et al.l 120091) and t he T6.5 brown 
dwar f 2MASS J10475385+2124234 (|Route fc Wolszczan 
20121 ). For simplicity we do not consider the 
dwarf binary since the behaviour of such a sys- 
tem could be significantly more complex than is 
represented by our simple axisymmetric model. 
Regarding the latter dwarf, two highly polarised 
bursts were observed at 4.75 GHz (implying a 
field strength of ^0.17 T), with spectral flux 
densities of 2.7 ± 0.2 mJy (where 1 Jansky = 
KT 26 W m~ 2 Hz" 1 ) and 1.3 ± 0.2 mJy, respec- 
tively. Assuming the radio emission is beamed 
into 1.6 sr and noting the object's 10.3 pc dis- 
tance, these flux densities imply a mean spectral 
luminosity of (3.3 ±0.7) MW Hz -1 . This observa- 
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Table 1: Properties of the UCDs considered in this paper. 



Property 


TVLM 513 


LSR J1835 


2M J0036 


Radius Rjjcd / Rj a 


1.005 


1.142 


0.927 


Rotation period t / h a 


1.96 


2.84 


3.08 


Polar ionospheric field strength Bi / T b 


0.3 


0.3 


0.17 


Distance d J pc a 


8.8 


5.7 


10.6 


Quiescent spectral luminosity L rq j MW Hz -1 c 


5.77 


2.77 


1.23 


Mean peak spectral flux density F r / mJy 


4.2 d 


2.4° 


0.5 f 


Pulse duty cycle At/t 


0.05 s 


0.1 c 


0.3 f 


Estimated peak spectral luminosity L r / MW Hz -1 h 


2.6 


0.8 


1.3 


Pedersen conductance Xp / mho 1 


0.5 


0.3 


0.5 



a Prom Table 1 of lHallinan et al.l J2008f ) 
kFrom lHallinan et al . (2008) 
c From Table 2 of lHallinan et al.l J2008fl 
d From Figure 1 of Hallinan ct al. | (120071') 
e From Figure 1 of lHallinan et al.l (120081 ) 
f From Figure 2 of lHallinan et al.l (120081 1 
9 From Figure 3 of lHallinan et al.l (120071 1 
''From equation 1151 
'Employed in Figure [4] 

tion is not included in Figure 2 since no rotation 
period is known, although it is likely to be at least 
2 h since only one burst was observed in each 
~2 h observing interval. Using the parameter 
values employed in this paper, our model indeed 
requires rotation periods of 2.1-2.8 h to produce 
radio powers consistent with the above luminosi- 
ties, such that further observations of this object 
are required to test this scenario. 

4. Summary and Discussion 

In this paper we have developed an axisym- 
metric model of the M-I coupling currents arising 
from departure from rigid corotation in the mag- 
netospheres of UCDs, based on that used previ- 
ously to study Jupiter's auroral oval. This sim- 
ple model generates radio emissions whose power, 
bandwidth, modulation period and duty cycle are 
consistent with those observed. The implication 
is that these UCDs are coupling with their space 
environments via magnetic fields, and that this 
coupling reduces the angular velocity of the high 
latitude ionopsheric regions, thus generating in- 
tense auroral emissions including radio emissions. 



In this study the angular velocity profile is 
used as the input to the model, whereas in real- 
ity the form of the profile will be dictated by the 
physical processes governing the magnetospheric 
flows. The magnetospheres of solar system plan- 
ets are driven by two sources of momentum, i.e. 
the solar wind and planetary rotation. Which 
of these is dominant depends on the individual 
circumstances, with Earth's magnetosphere being 
dominated by the solar wind, while Jupiter's and 
Satur n's are rotationally do minated. Previous 



work (jKuznetsov et al.l I2012T) has assumed that 



the radio from UCDs is emitted from field lines 
mapping to p e — 2Rucd, where p e represents the 
equatorial radial distance from the magnetic axis, 
on the basis that at this distance the gravitational 
energy density equals the corotational energy den- 
sity, such that it is hypothesised that beyond this 
distance any initially corotating plasma becomes 
centrifugally unstable and flo ws approximately ra- 
dially away from the object (|Ravi et al.ll201lf) . It 
is worth noting, however, that at Jupiter this dis- 
tance is also ~2 Rj, and yet the dominant plasma 
flow is azimuthal out to the magnetopause at ~40- 
100 Rj and no auroral emission is generated on 
field lines mapping to p e ~ 2 Rj . 
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A more complete equation of motion is 



Vp = j x B 



(16) 



where p is the plasma mass density, dv/dt is the 
acceleration of the plasma bulk flow with respect 
to inertial space, g is the acceleration due to grav- 
ity, p is the plasma pressure, j is the current 
density and B the magnetic field, and we note 
that in the corotating frame the dv/dt term com- 
prises the convective and local accelerations, cen- 
trifugal, Coriolis and differential rotation effects 
( Vasvliunasl 1983). The net result of the above 
force balance is that Jupiter's magnetosphere is 
radially distended into a magnetodisc structure, 
which doubles its size from that which might be 
expected from the pl anetary dipole alone (jCaudal 
19861 : Nichols 2011bl) . Plasma in Jupiter's mag- 
netosphere is, however, centrifugally unstable and 
diffuses r adially outward via flu x tube interchange 
motions (jBespalov et al.l 120061 ) before being lost 
down the tail via th e pinching off of plasmoids 
at around -100 Rj (|Vogt et al.l l2010h . If there 
were no ionospheric torque on the magnetospheric 
plasma its angular velocity would decrease with 
1/p 2 , where p is the distance from the spin axis, 
but the finite ionospheric Pedersen conductance 
allows a large-scale M-I coupling current system 
such as that discussed in this paper to flow, the 
j x B force of which balances the neutral drag in 
the ionosphere and maintains partial corotation 
throughout the closed magnetosphere. The equa- 
tion of motion which describes the angular velocity 
of the plasma under steady outflow fr om an inter - 
nal source was derived orig inally by EH! ljl97flh . 
and is given by 



Pe d 
2 dp e 



n 



UCD 



9. 



UCD 



AttT,' p F s \B^\ 
M 



\ ilucn j 



(17) 



where F e is the equatorial value of the the poloidal 
flux function F, related to the magnetic field B 
by B = (l/p)VF x 0, \B ze \ is the magnitude 
of the north-south magnetic field threading the 
equatorial plane, and M is the plasma mass out- 
flow rate. Analytic solutions to equation [T7] can 



be obtained, depending on the form of the mag- 
netic field, along with characteristic distances over 
which the plasma brea ks from rigid corotation 
(|Nichols fe Cowlevll2003l) . called the 'Hill distance' 
(not to be confused with the radius of the Hill 
sphere within which a body's gravitational force 
dominates). In the absence of knowledge of the 
magnetodisc structure of these UCDs, taking the 
dipole case represents a lowest order approxima- 
tion, for which the Hill distance is given by 



Ph \ 
Rucd J 



( 2nV* P B 2 R 2 



1/4 



q* L UCD 



M 



(18) 



where B eq is the equatorial field at the surface 
of the UCD. Employing £p = 0.5 mho and the 
canonical value of M — 1000 kg s _1 for the Io 
source rate, we then have pu ~ 775 Rucd- 



In comparison, the size of the magnetosphere 
Rmp is estimated by considering pressure balance 
between the magnetic field pressure of the com- 
pressed planetary dipole just inside the magne- 
topause and the magnetic, thermal, and dynamic 
pressures of the interstellar medium on the out- 
side, i.e. 



/ Rmp \ / '"m^eq 

\RUCD J \2po(pth + Pdyn + PB) 



k 2 Ft 2 
'-'e.q 



1/6 



(19) 

where k m is the factor by which the magnetopause 
field is enhanced by magnetopause currents, given 
by — 2.44 for a frontal boundary of realistic shape, 
i.e. lying between the values of 2 and 3 appropri- 
ate to planar and spherical boundaries, respec- 
tively (|Mead fe Beardlll964t \M exeevl 120051 ) . Tak- 
ing plasma number density, temperature and mag- 
netic field strength values for both the local inter- 
stellar cloud and the local bubble, along with a 
representative velocity for the UCD o f 50 km s" 1 
with respect to the local medium (|Vanhamaki 
201 lh , yields values of R mp ~ 713 and 820 Rucd 
for the local interstellar cloud and local bubble, 
respectively. Thus, on this basis corotation break- 
down may only be marginal inside the closed mag- 
netosphere, although a more detailed analysis of 
the force balance should be undertaken in the fu- 
ture to determine whether this would be the case 
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for a realistic current sheet field structure. 

It is also worth mentioning that while at least 
TVLM 513 may possib ly have a companion lo - 
cated at 20 - 50 R UCD (|Forbrich fc Bergerll2009l) . 
evidence for sate llite-induced radio emi ssion has 
not been found ( Kuznetsov et al.l 12012 ) , but in 



tosheath just upst ream of the magnetopause (e.g 



Siscoe et al. 1980). At planets in the solar system 



this regard we note that at Jupiter the bright- 
est and most significant auroral form resulting 
from the Io plasma source is the main oval, whose 
apparent modulation is at the planetary period, 
not Io's orbital period. If near-rigid corotation 
is then maintained throughout the closed magne- 
tosphere, an interesting point is that the linear 
velocities of the plasma near the magnetopause 
may reach a significant fraction of the speed of 
light, i.e. ~0.17c for rigid corotation with a pe- 
riod of 2 h at 800 Rucd- An order of magnitude 
approximation of the limiting distance at which 
plasma content can be maintained in corotation 
by inward magnetic stress is where the Alfven 
speed equals the corotation speed. The former is, 
of course, unknown since we do not have knowl- 
edge of the plasma density, but at high latitudes 
in Jupiter's magnetosphere it can approach the 
speed of light. If the limiting distance does lie sig- 
nificantly within the magnetopause then a super- 
Alfvenic radially-outward planetary wind may in- 
deed form, at which point the distinction between 
open an d closed field lines could become rather 
blurred (jKennel fc Coronitilll977l ). although it is 
worth noting that at Jupiter this distance is also 
within the magnetosphere and no swift radial out- 
flow occurs at this location. 

The amount of flux 'truly' open to the inter- 
stellar medium will depend on the reconnection 
rates on the leading and trailing sides of the flow 
with respect to the surrounding medium. Theories 
of magnetic reconnection have not yet reached the 
point where ab initio calculations of the rate of flux 
transport can be made, however it is known that 
reconnection occurs where a significant magnetic 
shear is present. Without in situ measurements 
we do not of course know what the interstellar 
magnetic field is at the magnetopause, but it is 
likely to be piled up against the boundary, since 
the ~50 km s" 1 proper motions typical of low 
mass stars are similar to the velocities observed in 
the shocked solar wind plasma in Jupiter's magne- 



the open-closed field line boundaries typically lie 
between ^10-15° co-latitude. The angular yeloc- 
ity of the open field line region ( Isbell et al.lll984h 
is given by 



n 



UCD 



Open 



1 + fl Y,* p V LIM 



(20) 



where vlim is the velocity of the UCD through the 
local interstellar medium, which for the above pa- 
rameter values yields (u> / Tin Cd) Open — 0.05, and 
may be much lower in the polar region if there is 
little aurorally-generated conductivity. Therefore, 
an open-closed field line boundary is a possible 
cause of the angular velocit y gradient discus sed in 
this paper. We note that ISchriiverl ( 20091 ) con- 
sidered the interaction of the magnetosphere with 
the interstellar medium, and ruled it out, by show- 
ing that the kinetic energy flux of the impinging 
medium is much less than typical X-ray and H-a 
luminosities. However, this does not consider the 
rotation of the body, which is a crucial property 
of these UCDs as discussed above. Here we pos- 
tulate that the interaction may instead be such 
that magnetic reconnection occurs on the front 
side of the magnetosphere and creates an essen- 
tially stagnant region of open field lines, whose 
angular velocities are given by equation 1201 and 
the transition from corotating closed flux to sub- 
rotating open field lines drives field-aligned cur- 
rents. This is the curre nt system described for 
Sat urn's auroras by e.g. ICowlev fc Buncel (|2008h 
and Cowlev et al.l ( 20041 ) and possibly a compo- 
nent of Jupiter's auroras poleward of the main 
oval (jCowlev et al.ll2003ll2005ll2007t ). Finally, it is 
worth noting that the auroral and radio emissions 
in our solar system caused by the above-discussed 
processes are all highly variable in terms of power, 
morphology and, in Saturn's case, modulation pe- 
riod, due to a number of reasons including vary- 
ing interplanetary medium parameters and inter- 
nal plasma source rates, such that it may be ex- 
pected that the radio emissions from UCDs are 
similarly variable. 
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